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A thermal equilibrium state of a quantum many-body system can be represented by a typical 
pure state, which we call a thermal pure quantum (TPQ) state. We construct the canonical TPQ 
state, which corresponds to the canonical ensemble of the conventional statistical mechanics. Only 
a single realization of the TPQ state gives all variables of statistical-mechanical interest. We show 
that TPQ states corresponding to different ensembles give identical thermodynamic results, and 
they are related by simple analytic transformations. We have thus completed the TPQ formulation 
of statistical mechanics. We also show that it has great advantages in practical applications. 
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Statistical mechanics is conventionally described by 
the ensemble formulation, in which an equilibrium state 
is represented by a mixed quantum state. Although 
its basic principles are introduced in the microcanonical 
ensemble, one can derive other ensembles, such as the 
canonical ensemble. All the ensembles are guaranteed to 
give identical thermodynamic results, i.e., they give the 
same equilibrium values of macroscopic variables. In this 
sense, all ensembles are equivalent. This makes statis- 
tical mechanics powerful and practical because for ap- 
plications the (grand) canonical ensemble is often more 
convenient than the microcanonical one. 

On the other hand, recent studies have shown that al- 
most all pure quantum states in a specified energy shell 
represent an equilibrium state [1-4]. Generalizing this 
fact, we have defined a thermal pure quantum (TPQ) 
state as a pure quantum state which represents an equi- 
librium state, and proposed a new formulation of sta- 
tistical mechanics based on the TPQ state [Hj. Since 
the TPQ state proposed there corresponds to the mi- 
crocanonical ensemble, we here call it a microcanonical 
TPQ state. Then, natural questions arise: Are there 
TPQ states which correspond to other ensembles? Are 
they equivalent? 

In this Letter, we construct a new type of TPQ 
state, the canonical TPQ state, which corresponds to the 
canonical ensemble. We show that it is equivalent to the 
microcanonical TPQ state, i.e., they give identical ther- 
modynamical results. We also show that they are related 
by simple analytic transformations. The TPQ states cor- 
responding to other ensembles, such as the grand canoni- 
cal ensemble, can also be constructed in a similar manner. 
These results complete the new formulation of statistical 
mechanics, according to which all quantities of statistical- 
mechanical interest can be obtained from a single real- 
ization of the TPQ state. This formulation is not only 
interesting as fundamental physics but also advantageous 
in practical applications because one needs only to con- 
struct a single pure state by just multiplying the hamil- 
tonian matrix to a random vector. The canonical TPQ 
state is particularly advantageous at low but finite tem- 



perature, whereas the microcanonical one is suitable for 
first-order phase transitions. As an illustration, we study 
the spin-1/2 kagome Heisenberg antiferromagnet (KHA) 
using the canonical TPQ state. 

Setup and definition - We use the same setup as that 
in Ref. That is, assuming an effective model that 
describes the energy scale of interest Q, we consider a 
discrete quantum system composed of N sites, which is 
described by a Hilbert space TLn of finite dimension A^, 
where A is a constant of 0(1) Q. [A = 2 for spin-i sys- 
tems.] To take the thermodynamic limit, we use quanti- 
ties per site, h = H /N and u = E/N, where H denotes 
the Hamiltonian and E the energy. We do not write ex- 
plicitly variables other than u and N , such as a magnetic 
field. We assume that the ensemble formulation gives the 
correct results, which are consistent with thermodynam- 
ics. 

Statistical mechanics treats macroscopic variables, 
which include 'mechanical' and 'genuine thermodynamic' 
ones. Mechanical variables, such as energy and spin- 
spin correlation functions, are the observables that are 
low-degree polynomials (i.e., their degree < to, where 
TO = o{N)) of local operators. In order to exclude foohsh 
operators (such as N^H), we assume that every mechan- 
ical variable A is normalized such that \\A\\ < KN™-, 
where K is a. constant independent of A and N. By 
contrast, genuine thermodynamic variables, such as tem- 
perature and thermodynamic functions, cannot be rep- 
resented as such operators. All genuine thermodynamic 
variables are functions of the entropy and its derivatives. 

We have defined the TPQ state for the case where 
it has a random variable(s), as follows 0]. A state \^) 

(e TLn) is called the TPQ state if {A)% ^ (A)] 



uni- 



formly as iV -H- oo, where {A)% = (iplAlip) / {^\^) , '^■' 



denotes convergence in probability, and (•)™^ the ensem- 
ble average. That is, for an arbitrary positive number e 
there exists a function r]^{N) that vanishes in the ther- 
modynamic limit and satisfies 



(1) 
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for every mechanical variable A. Here, P(a;) denotes the 
probability of event x. This means that for sufficiently 
large A'' just a single realization of gives the equilib- 
rium values of all mechanical variables. Note that |V') is 
never obtained by 'purification' of the density operator 
of the ensemble formulation, because {ip) is not a vector 
in an enlarged Hilbert space but a vector mUN- 

Canonical TPQ state - We have previously con- 
structed the microcanonical TPQ state, which is specified 
by {u,N) (i.e., independent variables are u, N) 15|. We 
now construct the canonical TPQ state, which is speci- 
fied by il3,N). [T = 1//3 is the temperature.] 

We take an arbitrary orthonormal basis of "Hjv, 
which can be a trivial one such as a set of product states. 
We also take random complex numbers {ci}i, which are 
drawn uniformly from the unit sphere X^i ^ 1- 
Then, = X^i'^iN) is a random vector in Hn- Note 
that this construction of |V'o) is independent of the choice 
of We will show that 



|/3,iV)=exp[-iVM/2]|V'o) 



(2) 



is the canonical TPQ state specified by {P,N). [It is 
unnormalized for later convenience.] We start with the 
rather obvious equality; 



,N\A\(3,N) _ J:^,JC*c,{^\e-^^^/^Ae-^^'^/^\J) 



W.N\(3,N) 
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where • • • denotes the random average, (A) 
TT[e-^^'^A]/Z{l3,N), and Z{I3,N) = Tre^^'^''. To in- 
vestigate how fast {A)f_N = ^1^1/5, ^)/(/?,A^|;3,iV) 
converges to this value, we calculate 



N 



[Af ^ ((i) 



TPQ _ / J\ can 



With the help of formulas in Ref . [8j , we find 
DN{Af < EN{Af, 

{{^Af)f-^ 



En{AY 



exp[7V/3{/(0;7V)-/(l//3;iV)}] 



(4) 

(5) 
(6) 



where ((Ai)^)-^^ = {{A - (A) f /(TjiV) is the 
free energy density, and we have dropped smaller-order 
terms. [We here use (T; N) instead of (T, N) to indicate 
that f{T;N) becomes independent of N in the thermo- 
dynamic limit.] Since {{AAf)^^^ < \\A\\^ < K^N^"", 
and NP{f{0;N) - /(l/^;iV)} = Q{N) at finite T (be- 
cause the entropy density s = —df/dT > 0), we find 
ENiAf < if^TV^^/e^^^), which becomes exponentially 
small with increasing N. We can also evaluate Dn{A)^ 
more strictly as 



DNiAf < 



213, N 



{(A) 



can 
2^, TV 



exp[27V/3{/(l/2/3;7V)-/(l//3;iV)}] 



(7) 



which is exponentially smaller than Ei\f{A)'^ for large N. 
On the other hand, a generalized Markov's inequality 
yields, for arbitrary e > 0, 



{AyrN >e)< D^iAf/e 



(8) 



The right-hand side < K'^N'^"' / €'^e^^^\ which vanishes 
exponentially fast with increasing N , for every mechani- 
cal variable A. Therefore, (^)J^ {AYp''% uniformly, 
which shows that |/3, N) is the canonical TPQ state. 

This result may be understood intuitively as follows. 
Although Ci 's distribute randomly, the law of large num- 
bers works well because an exponentially large number 
of Ci's are included in \P,N). As a resuh, A|/3, iV) 

and {/3^ N\/3, N) converge in probability to their random 
averages, whose ratio gives the correct result as Eq. ([3|). 
Although this intuitive discussion is easy to understand, 
it is important to obtain explicitly the upper bound of 
probability of error as formulas dS])-®. 

Genuine thermodynamic variables - We now show that 
1/3, N) also gives genuine thermodynamic variables cor- 
rectly. In the ensemble formulation, the partition func- 
tion Z{I3,N) gives / = -(1/;37V) In Z. Similarly, in our 
formulation the length of \P:N) gives / as 



(/3, A^|/3,iV) = Z(/3,iV)/A^ = exp(-7V/3/)/A 



(9) 



(|(/3,iV|/3,7V)/(/3,iV|AiV) 



> e 



< II exp[2A^/3{/(l/2/3; N) - /(1//3; A^)}]. (10) 

Therefore, (^,iV|/3,A^) Z(;3,iV)/A^, and a single re- 
alization of the canonical TPQ state gives /, with expo- 
nentially small probability of error, by 



1 

TV 



\n{(i,N\l3,N) -Kin A. 



(11) 



All genuine thermodynamic variables can be calculated 
from / [9|. Furthermore, using / obtained from this for- 
mula, one can estimate the upper bound of errors from 
formulas (H])-® and (|10l) . In this sense, our formulas are 
almost self-validating ones. This property is particularly 
useful in practical applications because one can confirm 
the validity of the results without comparing them with 
results of other methods . 

We can show similar results for the unnormalized mi- 
crocanonical TPQ state, 



\k) = {I - hf\i;o) (fc = 0,l,2,...), 



(12) 



which represents the equilibrium state at the energy den- 
sity Uk — {k\h\k) / {k\k) . Here, I is an arbitrary constant 
such that / > max{ eigenvalue oi h} Q. The squared 
leng th Qk = {k\k) gives the entropy density s{u;N) by 

m 

1 2k 

s{uk;N) ^ -lnQk-—\n{l-Uk)+lnX+0{l/N). (13) 



3 



This formula gives s more directly than the formula of 
Ref. f5l|, in which s is obtained by integrating /3. 

Equivalence of TPQ states and analytic relations- In 
the ensemble formulation, the equivalence of ensembles 
is crucial. Also in our formulation, the canonical and 
the microcanonical TPQ states are equivalent. In fact, 
Eqs. (dU and 1^ show that both TPQ states give 
the correct thermodynamic functions in the thermody- 
namic limit. Since /3f{l/f3) {— /3/(l//3; oo)) and s{u) 
(= s(m;oo)) are equivalent (because they are related by 
the Legendre transformation), so are both TPQ states. 

We can go one step further: We will show that these 
states are by themselves related by simple analytic trans- 
formations. To see this, we expand the exponential func- 
tion of e^'5'/2|/3,7V) = e^^('-'')/2|^o) as 



1/3, iV) 



fc=0 



(14) 



fe=0 



Here, j-^fe) = (1 A/Q/c)|^) is the normalized microcanon- 
ical TPQ state [3, and Ru = {N p/2fy/Q^/k\. We can 
prove that the coefficient Rk takes the maximum at k* 

(= e(iv)) such that [ni 



(5 = I3{uk*-N) + 0[1/N), 



(15) 



where /3(u; N) is the inverse temperature of the equi- 
librium state specified by {u;N), i.e., P{u;N) = 
ds{u;N)/du and itfc* = (V'fc* l^lV'fc*)- We can also prove 
that Rk vanishes exponentially fast for |fc — fc*| > Q{N). 
That is, Rk has a sharp peak at k* and relevant terms in 
the sum are localized in the range |fc — = o{N). This 
means that \I3,N) is almost composed of \ipk)^s whose 
temperature is close to 1//3. Moreover, we can prove that 
the sum is uniformly convergent on any finite interval of 
p. That is, if one fixes arbitrarily the upper bound /3max 
(> 0) of /3 depending on one's purposes, then the series 
is uniformly convergent for all /3 such that < /3 < /3max- 
We have thus shown that the canonical TPQ state can 
be expressed as superposition of the microcanonical ones. 
[Note that this expansion is obtained because we have 
expanded e^'^^'"''-'/^ rather than e~^^''/^.] Inversely, we 
can also obtain the microcanonical TPQ state from the 
canonical one, e.g., by 



1/3, /V) 



(16) 
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Different representations of the same equilibrium state 
- The canonical density operator e~^^^ /Z is invari- 
ant under time evolution. |/3,/V) is also invariant in 
the sense that it traverses various realizations of the 
same canonical TPQ state. To see this, let us inves- 
tigate e-^^'^^l^\fi,N) = exp[-iV/3/i/2]e-*^''*/'"'|i/'o). By 
taking the energy eigenstates {|7t)}n as the arbitrary 
basis used in the construction of |V'o), we have 



where = e~*^'^"*/''c„ and 
e„ is an eigenvalue of h. This is another realization of 
the random vector jt/'o)- Hence, e~^^^*'/^\f3^N) is another 
realization of the TPQ state for every t. Since almost all 
realizations of the TPQ state give identical results for 
macroscopic variables, as proved above, the TPQ state is 
macroscopically stationary in consistent with thermody- 
namics. This clearly illustrates the fact that the macro- 
scopic time-invariance of equilibrium states does not nec- 
essarily imply the microscopic time-invariance. More- 
over, according to experience, any quantum state rep- 
resenting an equilibrium state should be stable against 
weak external perturbations. We can show using the 
results of Refs. 12|, |l3| that the TPQ states do have 



such stability. These facts support that an equilibrium 
state can be represented by various microstates, includ- 
ing e'^P'' and \I3,N). 

Note also that one can generate in a quantum com- 
puter the TPQ state of any many-body system (as long 
as its properties are consistent with thermodynamics). 
This can be done, using Eqs. (|T2t and ([T4| (whose sum 
is terminated at fcterm; see below), only in a polynomial 
time Poly(/V). This observation may open new applica- 
tions of quantum computers. 

Practical formulas - It is practical to calculate |/3, A^) 
from |'i/'fe)'s through Eq. (1141) . because the microcanoni- 
cal TPQ states • ■ ■ , IV'fc) can be obtained iter- 
atively by simply multiplying (I — h) with I-^q) k times 
5]. Since Rk has a sharp peak at k* (given by Eq. (fTS)) ). 
one can terminate the sum at a finite number fcterm ■ It is 

sufficient to take fctorm such that fctorm — ^max = 0(-^), 

where fcj^ax is k* corresponding to /^max- Since we can 
show that k* = e{N) for any finite (3, k*^^^ = e(iV). 
Hence, fctorm = 6(/V). In this way, one can obtain \I3,N) 
by multiplying [l — h) repeatedly Q{N) times. 

All macroscopic variables can be calculated from the 
obtained \(3,N). One can also calculate them with- 
out obtaining \I3,N) explicitly. To see this, we note 
that all macroscopic variables can be obtained from 
{l3,N\A\l3,N) and (/3, Af|/3, /V), as shown by inequalities 
((HI) and (fTTI) . Since the latter is included in the former as 
the case oi A = 1, we consider the former. From Eq. (|14l) . 
{p,N\A\p,N) = e-^'5'Efc,fc'[W/2)'+''A!A:'!](fc|i|fc')- 
For the special case where [A, h] = 0, this reduces to 



{/3,N\A\(3,N) = {A}p^N, 



{A} 



I3,N 



E 



(2fc)! 

(A^/3)2fc+i 
(2fc+l)! 



(17) 

(fc|i|fc) 

(fc|i|fc + l). (18) 



Even when [A, h] 7^ 0, we can prove, rather surprisingly. 
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that Eq. ([nl) holds well. Specifically, 



N 



>e < 



En {Ay 



(19) 



This and inequality © imply {A}p^N/{^fi,N ^ (^))3'^w 
uniformly. Formula (1171) is useful because one needs only 
to calculate (fc|v4|fc) and (fc|yl|fc + 1) for all k < A:tcrm to 
obtain the results for all /3 < /3max- 

The TPQ formulation of statistical mechanics - It 
is straightforward to extend the above theory to the 
TPQ states corresponding to other ensembles, such as the 
grand canonical ensemble. We have thus completed the 
new formulation of statistical mechanics. It is summa- 
rized, for the microcanonical and canonical TPQ states, 
as follows. Depending on the choice of independent vari- 
ables, {E, N) or (/?, N), one can use either state, because 
they give identical thermodynamic results. A single real- 
ization of a TPQ state is sufficient for evaluating all quan- 
tities of statistical-mechanical interest. In addition, one 
can estimate the upper bounds of errors (which vanish 
as TV ^- oo) by formulas dS])-®, ^ and (HH). The mi- 
crocanonical and canonical TPQ states are transformed 
to each other by simple analytic relations, Eqs. (|T4| and 
P^ . Hence, getting either one implies getting both. Uti- 
lizing this fact, we have developed a practical formula 

Since the TPQ formulation is much different from the 
ensemble formulation, it will lead to deeper understand- 
ing of macroscopic quantum systems. It is also unique 
and advantageous as a numerical method: (i) It computes 
expectation values for a pure quantum state, whereas the 
other methods compute those for a mixed state (such as 
e~^ / Z). (ii) No limitation of models, i.e., applicable to 
any spatial dimensions and to complicated systems such 
as frustrated systems and fermion systems, (iii) Appli- 
cable to finite temperature. At finite T, there are an 
exponentially large number of states in the correspond- 
ing energy shell. This reduces accuracy of many other 
methods. By contrast, our method becomes more accu- 
rate as the number of relevant states increases, as explic- 
itly shown by formulas dS])-®, ^ and (HH). (iv) Al- 
most self-validating because one can estimate the upper 
bounds of errors from these formulas, (v) One can obtain 
the TPQ states by simply multiplying h repeatedly Q{N) 
times to a random vector. This is much faster, e.g., than 
the numerical exact diagonalization (ED) of the full spec- 
trum. [For example, it took only a few days to calculate 
all data in Fig. [T] on the Intel Xeon E5-2670 (2.6GHz).] 
(vi) To obtain the results for all /? < /3max, only two vec- 
tors, \k) and \k-\-l), are required to store in the computer 
memory, (vii) The orthogonality, {k\k') = for ^ A:', is 
not necessary at all. This is advantageous to large-scale 
computations. 

Regarding the choice between the canonical and mi- 
crocanonical TPQ states, one can use either depend- 
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FIG. 1. Specific heat c vs. T of the KHA for N = 18-30. 
(Lower Inset) / and s vs. T for = 30. (Upper Inset) The 
shapes of the clusters of A = 30 (left) and 27 (right). Those 
of 18a and 18b are shown in Ref. [l5| . 



ing on the purpose. For example, if one is interested 
in a first-order phase transition at which the specific 
heat c = du/dT = {dT/du)~^ diverges the microcanon- 
ical one is practically better, because T{u) is continu- 
ous (whereas u{T) is discontinuous) through the tran- 
sition. On the other hand, the canonical one is better 
when one studies low-temperature behavior of c, because 
du[T)/dT gets small {dT{u)/du diverges) as c — )■ 0. 

Application - As an illustration, we study the KHA, 
which is a frustrated two-dimensional quantum spin sys- 
tem. It was sug gested that c has double peaks at low 
temperature [IJ]. However, the problem is still in dispute 
due to the complexity of the frustration and the finite size 
effect [iM3|. We compute c, / and s for N = 18-30, tak- 
ing fcterm = 2000, for which the residual is evaluated to 
be less than 10"^"% for T > 0.02J. 

Figure [T] plots c which is obtained using d{h)p^^^/d(3 — 

-{Ch-{h)fj,?)p'j, ^ {A'},3,iv, where A = h-{h}p,N. 
We have also calculated c by using the difference method 
as d{h)f],/dp ~ i{h}f,+sp,N - {h}p^N)/Sf3. The differ- 
ence of these two numerical results is much smaller than 
the hue width of the data in Fig. [T] For A^ = 27 and 30, 
for which ED has never been performed, there is not a 
peak but a shoulder around T — 0.1 J, although the finite 
size effect may still be non- negligible. 

We also estimate the error caused from the random 
initial vector by using inequality ([S|). In the inset of 
Fig. [H we plot / (left scale), which are calculated from 
Eq. (ITT|) . Using the results for / and those for {(A^ — 
{A^}^^jv)^}/3,Af, we find that the (normalized) standard 
deviation Dm{A'^)/{A^}i3,n for A^ = 30 is less than 1% 
down to r — 0.1 J. The error of / itself is also estimated 
to be less than 1% down to T = O.IJ. Such a small error 
is attained because our method gets more accurate for 
larger entropy Ns and the KHA has relativel y la rge s at 
low temperature due to the frustration effect [16| . To see 
this quantitatively for A'^ = 30, we plot s = {u— f)l3 in the 
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inset of Fig. [T] (right scale). At T = 0.2 J there remains 
45% of the total entropy (= A'' In 2) . Such a large entropy 
makes Dj^{A'^) small. 

Finally, to confirm the validity, we compute c for N = 
18, for which the result of the ED is available 15|. For 
such a small cluster, the standard deviation estimated 
from inequality ^ is about 35% at T = O.IJ. Therefore, 
we have used Eq. ([3]) for = 18 only, taking average over 
100 realizations of the TPQ state. Our results (18a, 18b 
in Fig. [1]) agree well with those by the ED jl^ . 
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